% 频域相乘函数
function yn = myfilter(x,h)

N1 = length(x);
N2 = length(h);

N = N1 + N2 - 1;
XK = fft(x, N);
HK = fft(h, N);

YK = XK.*HK;
yn = ifft(YK, N);

if (all(imag(x)==0) && (all(imag(h))== 0))
    yn=real(yn);
end

end
